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1 Introduction 



At present, the investigation of W-pair production at LEP2 plays an important 
role in the verification of the Electroweak Standard Model (SM). Apart from the direct 
observation of the triple-gauge-boson couplings in e + e~ — > W + W~, the increasing 
accuracy in the W-pair-production cross-section and W-mass measurements has put 
this process into the row of SM precision tests [ [j]. The W-pair cross section is 
measured at the per-cent level, and the W-boson-mass determination aims at a final 
accuracy of 30 MeV. Experiments at a future e + e~ linear collider (LC) with higher 
luminosity and higher energy will even exceed this precision. 

To account for the high experimental accuracy on the theoretical side is a great 
challenge: the W bosons have to be treated as resonances in the full four-fermion 
processes e + e~ — > 4/, and radiative corrections need to be included. While sev- 
eral lowest-order predictions are based on the full set of Feynman diagrams, only 
very few calculations include radiative corrections beyond the level of universal radia- 
tive corrections (see Refs. [ 0, [| and references therein). These universal corrections 
comprise the leading process-independent effects; for e + e~ — > WW — > 4/ these in- 
clude universal renormalization effects (running or effective couplings), the Coulomb 
singularity at the W-pair-production threshold, and initial-state radiation (ISR) in 
leading-logarithmic approximation. The remaining corrections are usually viewed as 
non-universal and can only be included by an explicit diagrammatic calculation. In 
this article we describe the structure of the universal corrections in detail and discuss 
the size of the non-universal corrections for LEP2 und LC energies. This issue is not 
only theoretically interesting, it is also important in practice, since many Monte Carlo 
generators for W-pair production that are in use neglect non-universal electroweak 
corrections. 

The size of non-universal corrections was already estimated by inspecting the pair 
production of stable W bosons quite some time ago [ @, [§, |JJ- For LEP2 energies these 
effects reduce the total W-pair cross section at the level of 1-2%, but for energies in 
the TeV range the impact grows to 0(10%). For differential distributions the size of 
the non-universal corrections is usually much larger. In the following we investigate 
the corresponding corrections to off-shell W-pair production, e + e~ — > WW — > 4/, 
by inspecting total cross sections as well as angular and invariant-mass distributions 
with the Monte Carlo generator Racoon WW [ |J . 

2 Radiative corrections to off-shell W-pair production — 
state of the art 

Fortunately, to match the experimental precision for W-pair production a full 
one-loop calculation for the four-fermion processes is not needed for most purposes, 
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in particular for LEP2 physics. Instead it is sufficient to take into account only those 
radiative corrections that are enhanced by two resonant W bosons. For centre-of- 
mass (CM) energies Eqm not too close to the W-pair-production threshold, say for 
Ecm ~ 170 GeV, the neglected 0(a) corrections are of the order (a/ir)(Tw/Mw), 
i.e. below 0.5% even if possible enhancement factors are taken into account. The 
theoretically clean way to carry out this approximation is the expansion about the two 
resonance poles, which is called double-pole approximation (DPA). A full description 
of this strategy and of different variants used in the literature (some of them involving 
further approximations) can be found in Refs. [ @, 0, [§ fj ■ 

At present, two Monte Carlo programs include 0(a) corrections to e + e _ — > 
WW — > 4/ in DPA and further numerically important higher-order effects: 
YFSWW3 [ |7j and Racoon WW [ |j. The salient features of the two approaches, 
which are conceptually very different, as well as detailed comparisons of numerical 
results are summarized in Ref. [||. Further numerical results of the two programs 
can be found in Refs. [ ffDl O] . Both programs have reached an accuracy of roughly 



~ 0.5% for CM energies between 170 GeV and 500 GeV. For higher energies also 
leading electroweak two-loop effects become important (see also below). 

Figure |] shows a comparison of the results of Racoon WW and YFSWW3 with 
recent LEP2 data, as given by the LEP Electroweak Working Group [ [T^] for the 
Summer 2000 conferences. The data are in good agreement with the predictions of 
the two programs, which differ by about 0.3% at LEP2 energies. Below a CM energy 
of 170 GeV, the prediction in Figure [1| is continued by GENTLE [ [T|], which does 
not include the non-universal electroweak corrections. In its new version GENTLE 
is tuned to reproduce the DPA prediction of Racoon WW and YFSWW3 on the 
total cross section at LEP2 within a few per mill (see Ref. [ [3[]). 



3 Universal electroweak corrections — improved Born ap- 
proximation 

3.1 Preliminaries 

Universal radiative corrections are those parts of the full correction that are con- 
nected to specific subprocesses, such as collinear photon emission or running cou- 
plings, and lead to characteristic enhancement factors. Owing to their universality 
such corrections are often related to the lowest-order matrix element of the under- 
lying process. In the following we construct an improved Born approximation (IBA) 
for the processes e + e~ — > WW — > 4/ that is based on universal corrections only. For 
the production subprocess the IBA closely follows the approximation formulated in 
Ref. [ M for on-shell W-pair production. For the W decay the IBA is identical with 



the lowest-order prediction in the G„ scheme, as suggested in Ref. [ 14 . 
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Figure 1: Total WW production cross section at LEP2 as given by the LEPEWWG [ [T^] 



In order to define the IB A, we first need the lowest-order matrix element of the 



process 



e + (p+,o-+) +e (p_,<7- 



W + (/c + ,A + )+W-(A:_,A_) 

o-i ) + J 2 (k 2 , cr 2 ) + h(h, 03) + 7 4 (^4, <r 4 ) 



(3-1) 



The arguments label the momenta p±, ki and helicities <7j = ±1/2, Aj = 0, ±1 of 
the corresponding particles. The cross section that is defined by including only the 
so-called signal diagrams for W-pair-mediated four-fermion production, which are 
shown in Figure |2|, is called CC03 cross sectionP]. Note that the masses of the exter- 
nal fermions (not the ones in closed fermion loops) are neglected whenever possible. 
In the absence of photon radiation this, in particular, implies that we have helicity 
conservation for the initial e + e~ system, i.e. only the combination o~_ = — er + con- 



1 Of course, the CC03 cross section is a non gauge- invariant quantity. However, evaluated in the 
't Hooft-Feynman gauge it approximates the full cross sections for W-pair-mediated 4/ production 
very well, as long as no electrons or positrons are in the final state. Therefore, the CC03 cross 
section is widely used in the literature (see also Refs. [ 0, pi). 
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Figure 2: Lowest-order signal diagrams for e + e — ► WW — ► 4/ 



tributes, and we can define a = a_ = — cr+. For definite electron helicity a, the 
lowest-order CC03 matrix element is given by 



3 



M e Bo e -^^(p + , P ., k+, k-, k 2 + , k 2 _) = £ F; Born ( s , t)Mfa+,p-, k+, k-, k 2 + , kl 



n=l 



(3.2) 

where M a n are so-called standard matrix elements (SME) containing the spinor chains 
of the external fermions, and F° Born (s,t) are invariant functions containing couplings 
and propagator factors. In lowest order only three SME and invariant functions 
contribute. Following the notation and conventions of Ref . [ |6j , these read (u± = 
(l±7 5 )/2) 

M a 2 = v(p+)0 + - ^_)w ff «(p_)(e;e!.), 

Ml = v{p + )f + uj a u{P-W-K) - v{p+)^_u a u{p.){e* + k-) (3.3) 
with "effective W-polarization vectors" 



e 1 _ 



V2s w k\ — Myf + iM w r 



vv 



and 



e 2 



-^l,Born( S J — O 2+^ a 



J 

F 3 a Bom(s,t) = -F° Boin (s,t) = ?L - (3.5) 

The actual values of the input parameters e, Mw, M%, and s w depend on the input- 
parameter scheme. In the G^-scheme the electromagnetic coupling e is deduced from 
the Fermi constant using the tree- level relation e 2 = 4v^G>M{yS 2 tf , and the weak 
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mixing angle is fixed by the gauge-boson masses, which are independent input pa- 
rameters, = 1 — M^/Mf . 

Before we define the IBA we comment on the calculation of the full factorizable 
one-loop correction in DPAQ, which is described in Ref . [ [J , and its relation to the 
decomposition ( |3.2|) . In this case six independent SME contribute for each value of a, 
and the functions F£ contain standard loop integrals. Moreover, in order to guarantee 
the gauge invariance of the corrections, which is mandatory for consistency, it is 
necessary to perform an on-shell projection of the external fermion momenta fcj. This 
means that the ki are changed to related momenta hi in such a way that k 2 ± = Myj. 
The off-shell values k 2 ± are kept only in the propagator factors of (|3.4|) .P1 



3.2 Improved Born approximation 

The first step in the construction of the IBA consists in a modification of the Born 
matrix element in such a way that the universal renormalization effects induced by 
the running of a and by Ap are absorbed. This is achieved [ |J by the replacements 

e 2 

— -> 4v / 2G M M 2 ,, e 2 -> 4vra(s) (3.6) 



in the lowest-order functions F? Boin of (|3.5| ), which implies that weak-isospin exchange 
involves the coupling G^M^ and pure photon exchange the coupling a(s). The run- 
ning of the electromagnetic coupling is induced by light (massless) charged fermions 
only, i.e. we evaluate a(s) by 



a(M, 



«(*) = — (3-7) 



with the value a(M|) = 1/128.887 taken from the fit [ [HJ of the hadronic vacuum 
polarization to the empirical ratio R = a(e + e~ —>■ hadrons)/cr(e + e~ — ► fi + fi~). Thus, 
the basic matrix element for the IBA reads 

A<r^ WW ^ 4/ *>+,P-, k+, k-, k 2 + , k 2 _) = K,iba(s, t)Ml( P+ , P -, k+, k-, k 2 + , k 2 _) 

n=l 

(3.8) 



2 In DPA the virtual one-loop correction consists of factorizable and non- factorizable contri- 
butions. The factorizable corrections are the ones that are related to the W-pair-production and 
W-decay subprocesses. The non-factorizable corrections comprise the remaining doubly-resonant 
virtual corrections and include all diagrams with photon exchange between the production and 
decay subprocesses. 

3 This on-shell projection also renders the CC03 cross section gauge- invariant, leading to the so- 
called DPA Born cross section. However, the DPA Born cross section is a much worse approximation 
for the 4/ cross section than the CC03 variant (see also Refs. [ |L [1). 
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with 



*W, .) = - w. - -trtf «~ - (3 ' 9) 



Note that we have used the complex Z-boson mass Mf = Mf — iM w r z in order to 
regularize the Z resonance below the W-pair-production threshold; otherwise the ISR 
convolution over the reduced CM energy would lead to complications (see below). 

Another important virtual correction is induced by the Coulomb singularity near 
the W-pair-production threshold. We include this effect in the calculation of the 
"hard" IB A cross section <Tj B X ^ ww ^ 4 -^ 

/ d^r- ww ^( P+ ,p_) = 1 / d$ 4/ |A<r- ww - 4/ | 2 [l + 8 Coul (s, kl, kl) g{(3) , 

(3.10) 

where the correction factor <5coui is given by [ 0, [17 



5coui(s, fc + , fc_) = -^=r- Im <J In 



/3 ' I V/? + /? + A 







s 2 + fci + ki - 2skl - 2sk 2 _ - 2k\k 2 _ 



AM jj^m (3 . li; 

V s s 



with the fine-structure constant a(0). The auxiliary function 

g(P) = (l-f) 2 (3.12) 

restricts the impact of <5coui to the threshold region where it is valid. Its actual form 
(and its occurrence) is somewhat ad hoc but justified by a numerical comparison to the 
full 0(a) correction. Omitting this factor would lead to a constant positive correction 
of a few per mill above threshold, although the correct non-universal correction is even 
negative. 

The last ingredient in the IBA is the leading-logarithmic contribution induced by 
initial-state radiation (ISR). We follow the structure-function approach [ [TB|, where 
the full IBA cross section cxiba reads 

J dcr IB A = J dxi J dx 2 T^(x 1 ,Q 2 )T^(x 2 ,Q 2 ) J da^X ~* WW ~" 4f (x 1 p+,x 2 p-). 

(3.13) 
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The structure functions T^(x,Q 2 ) include the leading logarithms [a \n(Q 2 /m 2 )} n U P 
to order n = 3, and the soft-photon effects are exponentiated; the explicit expressions 
can also be found in Refs. [ ||, || . The QED splitting scale Q 2 is not fixed in leading- 
logarithmic approximation and has to be set to a typical momentum scale of the 
process. It can be used to adjust the IB A to the full correction, but also to estimate 
the intrinsic uncertainty of the IB A by choosing different values for Q 2 . 

Finally, we have to fix the W-boson width T-w in the evaluation of the IBA. In order 
to avoid any kind of mismatch with the decay, T w should be calculated in lowest order 
using the scheme. This choice guarantees that the "effective branching ratios", 
which result after integrating out the decay parts, add up to one when summing 
over all channels. Of course, if naive QCD corrections are taken into account by 
multiplying with (1 + a s /n) for each hadronically decaying W boson, these QCD 
factors also have to be included in the calculation of the total W width. 

Note that unlike the full one-loop calculation in DPA, the IBA is also applicable 
near the W-pair production threshold, since no pole expansion is involved. 



4 Comparison of the improved Born approximation with 
state-of-the-art results 

4.1 Total cross section 



In order to investigate the reliability of the IBA defined in ( 3.13 ), we have im- 
plemented this IBA in the Monte Carlo program RacoonWW, which provides 
state-of-the-art predictions for the full 0(a) corrections in DPA, as discussed above. 
For the following numerical evaluations we have adopted the input-parameter set of 
Refs. [1,1. 

Figure ||] compares different predictions for the total cross section (without any 
phase-space cuts) for the semileptonic process e + e~ — > udyU™I7 A1 (7) for CM energies 
-^cm up to 1 TeV. The IBA is evaluated for the two different scales Q 2 = s = 
E 2 CM and Q 2 = \t min \ = E 2 M (l - ft - 4M&/£g M ) /2 - M 2 ,, and "best" labels 
the RacoonWW prediction including all universal and non-universal corrections as 
described in detail in Ref. [ffl. The motivation for Q 2 = s is obvious; Q 2 = \t min \ is 
motivated by the fact that t m ; n corresponds to the minimal momentum transfer in the 
t-channel diagrams for forward scattering of on-shell W bosons, which dominates the 
cross section. The comparison of the corresponding relative corrections (normalized 
to the CC03 Born cross section in G M scheme) is shown in Figure |J 

For LEP2 energies, i.e. energies below 210 GeV, the difference between the two 
IBA versions reflects the typical uncertainty of 1-2% inherent in all predictions that 
neglect non-universal electroweak corrections. It turns out that the IBA with Q 2 = s 
is closer to the "best" prediction, with a maximal deviation at the upper LEP2 
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Figure 3: Predictions for the total cross section for the process e + e 
various approximations for radiative corrections 
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energies: ~ 0.6% at 200 GeV and ~ 0.8% at 210 GeV. Note that the "best" prediction 
is not included below 170 GeV, since the uncertainty of all predictions based on a 
DPA formally runs out of control near the W-pair-production threshold. On the 
other hand, the IB A does not suffer from this constraint. Since the IBA with Q 2 = s 
agrees with the "best" prediction near 170 GeV at the per-mill level this IBA version 
is an appropriate extrapolation of the "best" RacoonWW prediction down to the 
W-pair-production threshold. Of course, the theoretical uncertainty below 170 GeV 
is then of the order of one to a few per cent. 

For LC energies the IBA becomes more and more uncertain; for 1 TeV the two 
IBA versions differ already by ~ 5%. This signals that non-leading electroweak cor- 
rections become more and more important. The dominant effects are due to Sudakov 
logarithms [ [19|] of the type a In 2 (s/M^) which originate from the exchange of soft 
and collinear massive gauge bosons, i.e. W and Z bosons. The IBA does not account 
for these effects. Nevertheless the IBA with Q 2 = |i m i n | follows the "best" predic- 
tion within ~ 1-2% even for high energies. This is plausible, because the total cross 
section is strongly dominated by the t-channel pole for forward scattering for high 
energies, and this contribution is well approximated by the IBA. Note, however, that 
the good agreement could not be predicted without a comparison with the full DPA 
correction including non-universal electroweak corrections. On the other hand, it can 
be expected that the quality of the IBA with Q 2 = \t min \ also becomes worse if for- 
ward scattering is excluded or suppressed by phase-space cuts; this issue is further 
discussed below in the context of differential distributions. 



4.2 Differential distributions 

In order to define differential distributions, the kinematic information on the 
fermion momenta in e + e~ — > WW — > 4/(7) is required. In the presence of photon ra- 
diation, a consistent treatment of photons that are soft or collinear to charged fermions 
is crucial. If such photons are not recombined with the nearest charged fermion, i.e. 
if these photon-fermion systems are not treated as single "quasi-particles" , the bare 
fermion momenta in general lead to distributions that are not IR-safe, i.e. they in- 
volve mass-singular logarithms of the form alnrrif. For fermions other than muons 
such effects are definitely unphysical. In order to avoid such artifacts, we recombine 
soft and collinear photons according to the procedure^ described in Refs. [ ||, || . 

4 In this approach, first photons close to the beams are dropped in events, i.e. their momenta are 
set to zero. If the photon survives the cut to the beam, it is recombined with the charged fermion / 
if Mf 7 < M rec , where / is the fermion with the smallest invariant mass M| 7 = (pf + /c 7 ) 2 with the 
photon. Finally, events are discarded in which charged fermions are close to the beam. The size of 
the recombination cut M lec , thus, determines how many photons are recombined with the charged 
fermions. In Refs. [ |[ |) the two values M rec = 5 GeV and 25 GeV are chosen, defining a "bare" and 
a "calo(rimetric)" setup, respectively. 
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Figure 5: Predictions for the W + -production-angle distribution (left) and corresponding 
relative corrections (right) for the process e + e~ — ► ud^v^ at Eqm = 200 GeV based on 
various approximations for radiative corrections 
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Figure 6: Predictions for the W + -production-angle distribution (left) and corresponding 
relative corrections (right) for the process e + e~ — ► ud/x~F M at -Ecm = 500 GeV based on 
various approximations for radiative corrections 
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In Figures [| and || the full Racoon WW and IBA predictions for the W + -produc- 
tion-angle distribution are compared for the process e + e~ — > ud/x - ^ at the typical 
LEP2 energy of Eqm = 200 GeV and the LC energy Eqm = 500 GeV. The uncertainty 
of the IBA predictions induced by the QED splitting scale Q is about 1-2% and 
~ 5% for LEP and LC energies, respectively. The deviation of the IBA prediction 
from the full result is up to ~ 5% and ~ 5-10%, where the agreement is best for 
forward scattering (cos#\y+ — > 1), as anticipated above. The IBA uncertainty and the 
deviation from the full result further grow with increasing energy. This comparison 
is performed using the "calo" setup for the photon recombination of photons, but the 
sensitivity of the W-production-angle distribution to the recombination procedure is 
very weak (see Refs. [ [| 

The sensitivity to photon recombination is maximal in the invariant-mass distri- 
butions of the reconstructed W bosons, which can be seen by comparing Figures [?| 
and 0. The full correction shows a very strong dependence on the recombination pro- 



cedure, which was discussed in Ref. [ ITUJ in detail. The more inclusive recombination 
( "calo" ) leads to large positive corrections above resonance and thus to a shift of the 
resonance to the right, which can be of the order of some lOMeV [ [HI. Since this 



distortion of the W line shape is mainly induced by final-state radiation and radiation 
off the W bosons, the IBA, as defined above, does not account for this effect. It is 
obvious that the W-invariant-mass distributions can only be properly described if 
photon radiation from the W-decay processes is taken into account properly. Figures 
and H refer to the LEP energy Ecm = 200 GeV, but this conclusion is, of course, 
valid for all energies. 



5 Conclusions 

Electroweak radiative corrections to e + e~ — > WW — > 4/ typically amount to 
O(10%) at LEP2 energies and further increase for higher energies. We have ex- 
plicitly given analytical results for the universal process-independent corrections, 
which include effective coupling constants, the Coulomb singularity near the W-pair- 
production threshold, and leading ISR effects. They have been implemented in the 
Monte Carlo generator Racoon WW, which calculates the full 0(a) corrections in 
double-pole approximation. Using this program a comparison between universal ef- 
fects and the full correction has been presented. 

For LEP2 energies the universal corrections are dominant, and the remaining non- 
universal contributions reduce the total W-pair cross section by 1-2%. In angular 
distributions non-universal effects can reach several per cent, mainly in regions where 
the cross section is small. The radiative corrections to W-invariant-mass distributions 
lead to a distortion of the W resonance, which is mainly due to photon radiation off 
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Figure 7: Predictions for the W + -invariant-mass distribution (left) and corresponding rela- 
tive corrections (right) for the process e + e~ — ► ud/i^z/^ at Eqm = 200 GeV based on various 
approximations for radiative corrections, using the "calo" setup for photon recombination 
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Figure 8: Predictions for the W + -invariant-mass distribution (left) and corresponding rela- 
tive corrections (right) for the process e + e~ — > udfi~V^ at Eqm = 200 GeV based on various 
approximations for radiative corrections, using the "bare" setup for photon recombination 
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the charged final-state fermions and off the W bosons. This line-shape distortion is 
not accounted for by the above-mentioned universal effects. 

For LC energies, i.e. energies up to the TeV range, non-universal effects become 
more and more important. While the universal effects still describe W-pair production 
in the forward region within some per cent, non- universal corrections reach the order 
of several 10% for intermediate and large W-production angles. 
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